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The dynamics of the development of instability of the free surface of liquid helium, 
which is charged by electrons localized above it, is studied. It is shown that, if the 
charge completely screens the electric field above the surface and its magnitude is 
much larger then the instability threshold, the asymptotic behavior of the system can 
be described by the well-known 3D Laplacian growth equations. The integrability 
of these equations in 2D geometry makes it possible to described the evolution of 
the surface up to the formation of singularities, viz., cuspidal point at which the 
electric field strength, the velocity of the liquid, and the curvature of its surface 
assume infinitely large values. The exact solutions obtained for the problem of the 
electrocapillary wave profile at the boundary of liquid helium indicate the tendency 
to a charge in the surface topology as a result of formation of charged bubbles. 
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I. INTRODUCTION 

It is well known that the liquid helium surface may be charged to high values of 

the surface density of a negative electric charge. This is due to the fact that, on the one 
hand, electrons are attracted to the surface by weak electrostatic image forces and, on the 
other hand, the liquid helium boundary is a potential barrier for electrons, which prevents 
their penetration in the bulk. An important feature of liquid helium as a dielectric with a 
low polarizability is the relative weakness of the image forces, as a result of which the mean 
distance between localized electrons and the surface is much larger than the atomic spacing. 
Consequently, the electrons are not bound to individual atoms of the substance and form a 
two-dimensional conducting system. 

The ability of electrons to move freely over the surface of liquid helium ensures the 
equipotential nature of this surface over characteristic hydrodynamic times and scales. A 
charged surface of a conducting liquid also possesses this property, the only difference being 
that the electric field cannot penetrate into a conducting medium, while liquid helium is 
not subjected to such a limitation. This enabled Gor'kov and Chernikova to extend 

a number of classical results from the theory of instability of a liquid metal surface in 
an external electric field (5, |fi, [?| to the case of the charged boundary of liquid helium 
(the geometry of the system is shown schematically in Fig. |IJ). For example, a natural 
generalization of the dispersion relation for linear waves on the surface of a conducting 
liquid is the following dispersion relation for liquid helium: 



u 2 = gk+-k 3 - ^ — k\ (1) 
p 4vrp 

where u is the frequency, k is the wave number, g is the acceleration due to gravity, a is the 
surface tension, p is the density of the medium, and E' and E are the electric field strengths 
above the liquid and in the bulk of it, respectively (E = for a conducting medium). It 
follows hence that for 

E' 2 + E 2 < E 2 = Sn^gap 

the inequality uo 2 > holds for any k and, hence, small perturbations of the surface do not 
build up with time. In the case when the sum of the squares of the fields E' 2 + E 2 , which 
plays the role of an extrinsic controlling parameter, exceeds the critical value E c 2 , a region 
of wave numbers k for which uj 2 < is formed. This corresponds to an aperiodic instability 
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FIG. 1: Schematic diagram of the surface of liquid helium, charged by electrons, in a parallel-plate 
capacitor. 

of the liquid boundary. 

The buildup of perturbations of the surface inevitably transforms the system to a state 
in which its evolution is determined by nonlinear processes. The nature of their effect 
can be estimated most easily in the vicinity of the instability threshold, i.e., for a small 
supercriticality e = (E 2 + E' 2 — E c 2 )/E c 2 , when only perturbations with wave numbers 
close to ko = y/gp/a increase and we can pass to envelopes in the equations of motion. 
For example, Gor'kov and Chernikova |8j proved that, in the case of 2D symmetry of the 
problem, the complex amplitude A(x,t) of perturbation of the surface obeys the nonlinear 
Klein-Gordon equation 



(gko)- 1 A tt = 2eA + kf A xx + (2S 2 - 5/8) A\A\' 



(2) 



where S = (E 2 — E' 2 )/ E 2 is the dimensionless parameter characterizing the surface charge 
density. It can be seen from Eq. (2) that, depending on the value of parameter S, the 
nonlinearity either saturates the instability, or, conversely, facilitates a burst of the pertur- 
bation amplitude. A similar conclusion can also be drawn in the general (3D) case with a 
correction taking into account the fact that the nonlinearity in the first nonvanishing order 
plays a destabilizing role due to the interaction of three waves forming the hexagonal struc- 
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lire. As in 2D case, cubic nonlinearities produce a stabilizing effect for small values of S 
^, [it]]. Consequently, for a low surface charge density (when the values of E and E' are 
close), a steady-state relief of the liquid helium boundary may be formed. In this case, the 
standard perturbation theory in the small parameter, viz., the characteristic slope of the 
surface, can be used for studying the structures being generated (see [l]] and the literature 
cited therein). 

The processes occurring in the supercritcal region of electric fields and for relatively 
large electron surface charge screening the field above the liquid surface to a considerable 
extent have not been investigated in detail theoretically. This is due to the fact that, 
in these cases, the development of instability violates the small-angle approximation. For 
example, the analysis of the behavior of the charged boundary of liquid helium by high-speed 
microphotography carried out by Vololdin et al. [12;] proved that the dimples appearing on 
the surface are sharpened over a finite time (the bubbles which are subsequently formed at 
the tips carry the charge from the helium surface to the positive plate of the capacitor). 
In view of the considerable nonlinearity of such processes, their description requires the 
construction of solutions to the fundamental equations of the electrohydrodynamics of liquid 
helium. 

In the present work, it will be shown that, when the condition E 3> E', which corresponds 
to complete screening of the field above the liquid by the surface electron charge, is satisfied 
along with the condition that the electric field strength considerably exceeds the critical 
value, E 3> E c , the equations of motion of liquid helium have an infinitely large number of 
exact analytic solutions. Their analysis has facilitated a considerable advance in the analysis 
of unsolved problems in the electrohydrodynamics of liquids with a free surface, which are 
associated with the formation of singularities (cusps) and with considerable changes in the 
surface geometry (formation of bubbles) . 

In Section 2, the equations of a vortex-free flow of liquid helium with a free surface 
charge are considered. In the limit of a strong electric field, when the effect of the force 
of gravity and capillary forces can be neglected, the approach to an analysis of the liquid 
helium dynamics proposed in our earlier work |l3( is developed. This approach is based on 
the separation of two branches corresponding to solutions increasing and decreasing with 
time in the equations of motion. In Section 3, it is shown that the asymptotic behavior of 
the system is given by the well-known equations describing the Laplacian growth in the 3D 
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geometry (the motion equipotential boundary with the velocity determined by the normal 
derivative of the harmonic potential). Section 4 is devoted to an analysis of the dynamics of 
the formation of cuspidal dimples on the helium surface in 2D geometry, when the Laplacian 
growth equations have an unlimited number of exact nontrivial solutions. The propagation 
of nonlinear surface waves in the short-wave region in which the surface pressure must be 
taken into account along with the electrostatic pressure is considered in Section 5. It is 
shown that the problem of the profile of a progressive electrocapillary wave at the liquid 
helium boundary has exact analytic solutions similar to the Crapper solutions for capillary 
waves ^| . These solution are used for obtaining a nonlinear dispersion relation for surface 
waves of an arbitrary amplitude, whose analysis led to a number of conclusions concerning 
the stability of the charged surface of liquid helium to finite-amplitude perturbations and 
the domain of the existence of wave solutions to the electrohydrodynamic equations. In 
Section 6, the simplest axisymmetric solutions of the equations of motion, describing the 
pulling of the surface into the bulk of the liquid at a constant rate, are analyzed. 



II. INITIAL EQUATIONS: THE LIMIT OF A STRONG FIELD 

Let us consider the potential motion of an ideal dielectric liquid (liquid helium) with a 
free surface charged by electrons in an electric field. We assume that, in the unperturbed 
state, the boundary of the liquid is a flat horizontal surface z = and the field vector is 
directed along the z axis of our system of coordinates (Fig. Q). We introduce a function 
r](x,y,t) specifying the deviation of the boundary from the plane. Then, the shape of the 
perturbed surface of liquid helium is described by the equation z = i](x,y,t). The velocity 
potential <3> for an incompressible liquid satisfied the Laplace equation 

V 2 $ = 0, (3) 

which must be supplemented with the dynamic boundary condition 

^t + L V L = ^ J1L + -Vx- % = -97], z = r)(x,y,t), 4 

where <p is the electric potential in the liquid (we assume that the charge completely screens 
the field above the helium surface). The first term on the right-hand side of the time- 
dependent Bernoulli equation (4) is responsible for electrostatic pressure, the second is re- 
sponsible for capillary pressure, and the third takes into account the effect of the field of 
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gravity. We assume that the characteristic spatial of surface perturbations is smaller than 
the size of the region occupied by the liquid. In this case, we can write 

$ — > 0, z — > -oo, (5) 

i.e., the motion of the liquid attenuates at infinity. The time evolution of the free surface is 
determined by the kinematic relation (the condition that the liquid does not flow through 
its boundary): 

rjt = <$> z - V±r] ■ Vl$, z = r)(x, y, t). (6) 
Finally, the electric potential <p in the absence of space charges satisfies the Laplace equation 

W = 0, (7) 

which must be solved under the condition that the liquid helium boundary is equipotential 
and the field is uniform at an infinitely large distance from the surface: 

(p = Q, z = r](x,y,t), (8) 



ip — > —Ez, z — > — oo. (9) 

It should be noted that, in zero electric field (E — and, hence, Vy? = 0) the above equations 
coincide with the equation of motion for a thick layer of liquid in the field gravity. 

Let the electric field strength exceed considerably its critical value (E ^> E c ), and let the 
following relation hold for the characteristic wavelength A of surface waves: aE~ 2 C A < 
E 2 /(gp). It follows from the dispersion relation (1) that, in an analysis of small-amplitude 
surface perturbation, we can disregard the effect of both the capillary forces and the force of 
gravity. In Secton 4, we will prove that this statement also holds for finite-amplitude surface 
perturbations. This means that we can omit the last two terms on the right-hand side of 
the boundary condition (4) and take into account the electrostatic pressure alone. 

Now we pass to the dimensionless notation, assuming that the unit of length is equal to 
A, the unit of electric field strength is E, and the unit of time is A-E _1 (47rp) 1//2 . In this case, 
the equations of motion (3)-(9) assume the form 

vV = 0, V 2 $ = 0, (10) 
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$ t + (V$) 2 /2 + (V^) 2 /2 = 1/2, z = rj(x, y, t), 



T) t = $ z - V±T] • Vj.$, 2 = 77(x, y, t), 



(p = 0, z = r)(x,y,t), 



$ -> 0, 



z — > — oo, 



(11) 

(12) 
(13) 
(14) 



V? — > — z, z — > — oo. (15) 

Let us write these equations in the form which does not contain function r\ explicitly and 
introduce the perturbed harmonic potential Cp — ip + z attenuating at infinity (ip — > as 
z — > — oo). At the boundary, we have <p\ z=T) = rj. This readily leads to relations 



Vt 



l-<Pz 



Z=7) 



z= V 1 - ^Z 

which allow us to eliminate f] from Eq. (12). The kinematic and dynamic boundary condi- 
tions (11) and (12) can be transformed to 

&-$* = -V£-V$, z = r)(x,y,t), 

*t-<Pz = -(V$) 2 /2 - (Vy3) 2 /2, z = V (x, y, t). 
Adding and subtracting these equations, we obtain 

(<p + $) t -{<p + $) 2 = -(V(£ + $ )) 2 /2, * = 77(x, y, t), 

(<p - $) t + (<p- $ ) z = +(V(v3 - $)) 2 /2, z = r)(x, y, t), 

i.e., the boundary conditions can be specified separately for the sum and the difference of 
the harmonic potentials <p and $. It is convenient to introduce a pair of auxiliary potentials 

(f>W(x,y,z,t) = (<p±$)/2. 

Using these potentials, we can write the equations of motion in the following symmetric 
form: 

V 2 (±) = 0, (16) 
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4 ±} = ±4>F> T (V0W) 2 , z = V (x, y, t), (17) 

(±) -> 0, ^ -> -oo, (18) 
while the shape of the liquid helium boundary is determined from the relation 

V= (0 (+) + ( " ) )| 2= ,- (19) 

Thus, the equations of motion can be split into two systems of equations for potentials 
<// + ) and 4>^~\ the relation between which is given by the implicit equation for the shape 
of the surface (19). It is important that these equations are compatible with the condition 
0*") = or with the condition = 0. In the next section, we will show that the former 
condition corresponds to the solutions of the problem whose amplitude increases with time, 
while the latter (which is of no interest to us), to damped solutions. 

The possibility of separating equations into individual branches is due to the symmetries 
of the elect rohydro dynamic equations, which can be easily seen when the Hamilton formalism 
is used. Indeed, the equations of motion (10)-(15) for a liquid with a free surface possess a 
Hamilton structure, the function r](x, y, t) and ip(x, y, t) = $| 2=r? being canonically conjugate 

3, 



quantities 



5H SH 

or) dip 



where the Hamiltonian H coincides to within constants with the total energy of the system: 

H = K + P, K = J d 3 r, 

z<r) 

P = J 1 -Wf T = -JSi£*r. 

z<r] z<r) 

It should be recalled that the harmonic potentials $ and <p attenuate for z — > — oo and 
their values on the surface are defined by the functions ip and rj, respectively. Consequently, 
if ip — t], then $ = and the kinetic energy functional K coincides, expect for the sing, 
with the potential energy functional P. This allows us to write the Hamilton equations of 
motion using the functional K alone: 
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It can be seen that, if we set ip = 77 in these equations, they will coincide. This means that 
the condition ip = 77 or (which is the same) the condition 0^ = is compatible with the 
equations of motion for liquid helium. Similarly, we can prove that the Hamilton equations 
coincide for tp = —i], which corresponds to the condition = 0. It should also be noted 
that the equations describing the evolution of the system on the branches = and 
<p(~} = coincide except for the substitution t — > —t, which is associated with the time 
reversibility in the Hamilton equations of motion. In this case, the conditions <f>^ = 
single out the solutions of the problem for which H is equal to zero. 

III. INCREASING BRANCH: STABILITY 

In the linear approximation whose applicability is limited by the condition of the smallness 
of the slopes of the surface |Vj_?7| <C 1, the boundary conditions (17) assume the form 

4 ±) = ±<f>^\ z = 0, 

and Eqs. (16)— (19) split into two independent systems. The dispersion relations for these 
systems can be found by substituting potentials in the form ~ e kz+ikr ± -iu>t _ ^his gives 

= ±ik 

(the same result follows directly from the dispersion relation (1) considered in the strong 
field limit). It can be seen that, for one branch, small periodic perturbations of the surface 
increase exponentially with the characteristic times while, for the other branch, these 
perturbations attenuate. In this case, for large periods of time, we can assume that = 
and consider only equations for potential (p( + \ Let us prove that this statement is also valid 
in the general case, when the evolution of the surface is described by nonlinear equations 
(16)-(19). 

We assume that, in the nonlinear equations of motion (16)-(19), 

0(+) = ip + z, = 0, 

which, in accordance with the results of linear analysis, isolates the solutions increasing with 
time. Passing to the moving fame of reference {x, y, z'} = {x, y,z — t} in which the plane 
unperturbed surface of the liquid moves downwards (i.e., in the direction opposite to the z' 
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axis) at a constant velocity, after simple transformations, we obtain 



VV = 0, (20) 

r/ t = d n <p Vl + (V±77') 2 , z' = rf(x,y,t). (21) 

^ = 0, z' = r)'(x,y,t) (22) 

V? — > — z', z' — > — oo, (23) 



where ^'(x, y,t) — rj — t and <9 n denotes the derivative along the normal to the boundary of 
the liquid. These equations define explicitly the motion of the free charged surface of liquid 
helium z' = r)'(x,y,t). They coincide with the equations describing the so-called Laplacian 
growth, viz., the motion of the phase boundary with a velocity directly proportional to 
the normal derivative of a certain harmonic scalar field (ip in our case). Depending on the 
chosen frame of reference, this field may have the meaning of temperature (Stefan's problem 
in the quasi-stationary limit), electrostatic potential (electrolytic deposition), or pressure 
(flow through a porous medium). 

Let us prove that the solutions of Eqs. (10)— (15) corresponding to system (20)-(23) are 
stable to small perturbations of potential cf)^~\ It should be noted that the motion of the 
liquid boundary described by Eqs. (20)-(23) is always directed inwardly; this is associated 
with the principle of the extremum for harmonic functions. Let function rj' at the initial 
instant t = be a single-valued function of variables x and y. In this case, for t > 0, the 
following inequality holds: r)'(x,y,t) < r)'(x,y,0). In the original notation, we have 

V(x,y,t) <r](x,y,0) + t (24) 

for any x and y. This inequality remains valid for small perturbations of 0^ also, when 
the effect of potential in relation (19) can be disregarded as compared to the effect of 
potential 4>^ + \ and the motion of the boundary is described by the same Eqs. (20)-(23). 

As regards the evolution of potential 4>^~\ it is described, for small |V0^ _ ^|, by Eqs. (16)- 
(18), where it is sufficient to consider the condition (17) at the boundary in the linear 
approximation: 

<PV = z = rj(x,y,t). 
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Let us suppose that, at the initial instant t — 0, the potential distribution is described by 
the following expression: 

(_) |t=o = </>o(x,y,z), 

where 0o is a certain function which is harmonic for z < r)(x,y,0) and attenuating for 
z — > — oo. In this case, the temporal dynamics of potential is described by the expression 

0(-) = (f) ( x ,y,z- t). 

It can be seen from this expression that the singularities of the function are displaced 
in the direction of the z axis and can exist only in the region 

z> V {x,y,0) + t. (25) 

A comparison of this inequality with (24) shows that the singularities of potential (f)^ do 
not approach the liquid helium boundary z = i](x, y, t) and, hence, the value of the potential 
at the surface does not increase with time. It should be noted that, otherwise, the solutions 
obtained by us for would be inapplicable. 

In view of incompressibility of the liquid, the level of its surface (the value of function rj 
averaged over the spatial variables) is not displaced. On the other hand, the boundary of 
the region defined by inequality (25) and averaged over x and y, in which singularities of the 
function <j>^ occurs, moves upwards at a constant velocity. This means that the singularity 
moves away from the surface of liquid helium and the perturbation of 0*") relaxes to zero. 

Thus, we have proved that, as t — > oo, we have 

<p(x,y,z,t) + z-> $(x,y,z,t), 

and Eqs. (20)-(23) describe the asymptotic behavior of liquid helium with a charged surface 
in a strong electric field. 

IV. SOLUTIONS OF 2D EQUATIONS OF MOTION 

In the previous section, we proved that the analysis of the 3D potential motion of liquid 
helium in a strong electric field can be reduced to analysis of Eqs. (20)-(23) describing the 
three-dimensional Laplacian growth. The exact solvability of these equations in the 2D 
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geometry will allow us to effectively study the dynamics of the development of instability of 
the charged surface of a liquid, including the formation of singularities in it. 

We assume that, in the system of equations (20)-(23), all quantities are independent of 
variable y (variable y'). We introduce the function w — v — itp of the complex argument 
Z = x + iz', which is analytic for z' < r]'(x, t) (this is the co-called complex potential of the 
field correct to a constant factor). Here, v is a function harmonically conjugate to ip and 
such that the condition v = const defines the electric field lines in the medium. Clearly, 
w^ZasZ-^x — ioo. 

It is convenient for the subsequent analysis to pass to a system of coordinates in which 
the role of the independent variable is played by quantity w and the role of the unknown 
function is played by function Z which is analytic in the lower half-plane of the complex 
variable w (i.e., for if > 0). It follows from condition (23) that the following condition holds 
at infinity: 

Z — > w, w — > v — ioo. (26) 

We can also obtain the condition for Z at the boundary if — of the half-plane. The profile 
of the liquid helium surface can be specified by the parametric relations 

z' = z'(v, t) = r)'(x(v, t),t), x = x(v, t). 

Using these relations, we can easily express the normal velocity of the surface and the electric 
field strength appearing in formulas (21) in terms of the functions z'(v,t) and x(v,t): 

Vt _ z t x v ~ x t z v q 1 

Substituting these relations into the condition (21) at the surface, we obtain 

Z t X v XfZ v 1, 

or, which is the same, 

lm(Z*Z w ) = 1, w = v. (27) 

Thus, we arrive at the problem of determining the function Z, which is analytic in the 
lower half-plane of the complex variable w and satisfies conditions (26) and (27). The 
nonlinear condition (27) is the so-called Laplacian growth equation which is widely used 
for describing the 2D motion of the boundary between two liquids with noticeably different 
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viscosities 



17( , the evolution of the free surface of a liquid in the field of gravity |l8|, 3| , 



and so on. The Laplacian growth equation is int egra ble in the sense that it has an infinitely 
large number of particular solutions of the form |20J] 

N / N \ 

Z{w) = w — it — i a n In (w — w n (t)) + i I a n J In (w — wo(t)) . (28) 

n=l \n=l / 

Here, a n are complex constants, and the functions of time w n satisfy the condition lm(w n ) > 
(singularities of the function Z can only be in the upper half-plane of of the complex 
variable w). The last term in expression (28) was supplemented to ensure the fulfillment of 
condition (26) and, hence, the condition of localization of the perturbation of the surface in 
a certain region: rj — > for \x\ — > oo. We can set Im(wo) ^ im (^n); i n this case, the effect 
of this term on the evolution of the surface is negligibly small. 

Substituting expression (28) into Eq. (27) and decomposing the obtained expression into 
simple fractions, we obtain a system of N ordinary differential equations for w n (t): 

N . . „ 

< — r = °- 

Integration with respect to t leads to the following iV transcendental equations: 

TV 

W n + it + i ^ a m ln ( W n ~ W* m ) = C n , 
m=l 



where C n are arbitrary complex constants. 

Let us consider the simplest solutions of this type, which correspond to N — 1, Re(iUi) = 
and a\ = ±1: 

Z(w)=w-it^ihi(w-iq(t)), (29) 

q(t)±]nq(t) = l + t c -t, (30) 

where q = Im(wi) and t c is a real constant. The form of a solitary perturbation corresponding 
to Eqs. (29) and (30) is specified by the parametric expressions 

z(v, t) = z'(v, t) + 1 = =p ln ^/v 2 + q 2 (t), 



x(v, t) — v ± arctan (v/q(t)) 
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FIG. 2: The profile of the liquid helium surface, corresponding to the "one-finger" solution of the 
Laplacian growth equation; a\ = 1, q = 10 -4 , and u>o = «5. 

Let us suppose that a% = +1 and we are dealing with a solitary perturbation of the 
surface, which is directed "upwards". It can be seen from Eq. (30) that for large values 
of t, the quantity q ~ e~* and, hence, the surface perturbation amplitude increases linearly 
with time: z\ v= q — > t as t — > oo. This is the "one-finger" solution of the Laplacian growth 
equation (see Fig. |2J). It can easily be proved that similar solutions are possible in the 3D 
case also. It can be seen from Eqs. (20)-(23) describing the three-dimensional Laplacian 
growth that, if the surface initially contains a region in which the field strength d n ip is small 
(e.g., in the vicinity of the apex of a 3D fingerlike perturbation of the surface), its velocity in 
the coordinates {x, y, z'} is also small. In the laboratory reference frame, this corresponds 
to a jet flowing at a constant velocity in the direction of the z axis. 



15 




FIG. 3: The profile of the liquid helium surface at the instant of formation of a singularity (cusp); 
ctl = —1, q = 0.8, and wq = «4. 

Let us now consider a solitary perturbation of the surface, which is directed "downwards" 
(ai = —1, q(t) > 1). This solution exists only during a finite period of time, leading to the 
formation of a singularity on the liquid surface, viz., cuspidal point of the first kind (Fig. OJ), 
at instant t = t c . Indeed, expanding z and x into power series in v and r taking into account 
the fact that the function q(t) in the vicinity of the t c satisfies the relation 



?(*) 



2r, 



T = t c — t, 



we obtain the following expressions in the main order: 

z = v 2 /2 + V2t, x = v 3 /3 + vV2t. 



(31) 



It can be seen that, at instant r = (i.e., for t = t c ), the shape of the surface in the vici nity 
of a singularity point is defined by the relation 2z = \3x\ 2 ^ 3 , which corresponds to a cusp [31]. 
It was indicated in 17, 23] that the singularities of z 3 ~ x 2 are general-position singularities 
for processes described by the Laplacian growth equation. Similar solutions of the equations 
of motion for liquid helium with the charged boundary reflect the experimentally observed 
tendency 0,0) to the emergence of dimples on the surface, which become sharpened over a 
finite time. From the mathematical point of view, the emergence of singularity on the liquid 
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surface is associated with vanishing of the Jacobian of the transformation {x, z'} — > {v, ip} 
for ip — v — t — 0. At a cusp, the electric field strength increases indefinitely along with 
the velocity of the surface over a finite time interval: 

It is important to note that the singular solution of the problem described by expressions 
(31) is also valid in the case when the field above the surface is not screened completely; i.e., 
the condition E' <C E does not hold. As a matter of fact, in the vicinity of a singularity, 
the condition of the smallness of the field above the surface as compared to the field in the 
bulk of the liquid naturally holds. In addition, the condition A <C E 2 / (gp) is not necessary. 
This is due to the fact that the amplitude of surface perturbations remains finite, and the 
effect of the gravity forces is always negligibly small in the vicinity of the cusp. 

Let us now consider the capillary effects. The surface and the electrostatic pressure in 
the vicinity of a singularity can be estimated easily: 

P,s ~ ar] xx ~ ap 1/2 E- X T~ x , P E ~ (Vy?) 2 ~ Ap 1/2 Sr _1 . 

Here, we have returned to the dimensional notation. It can be seen from these expression 
that, when the condition A aE~ 2 is satisfied, the capillary forces are small as compared 
to electrostatic forces and, hence, can be disregarded at the stage of formation of cusps. 
This is the only necessary condition of the applicability of the Laplacian growth equation 
and its solutions (31) in the vicinity of singularities. 

V. ELECTROCAPILLARY WAVES 

Let us consider the case when the characteristic length of the surface waves is comparable 
with the value of aE~ 2 and the capillary effects must be taken into consideration. We 
assume that condition E 3> E c is satisfied; in this case, the effect of the force of gravity can 
be neglected. The dispersion relation (1) for electrocapillary waves at the charged boundary 
of liquid helium for E' = in the dimensionless notation introduced in Secton 2 assumes 
the form 

uj 2 (k) = k 3 - k 2 , (32) 

where the value of A = AnaE~ 2 is taken for unit length. It can be seen from Eq. (32) 
that uj 2 < for k < 1 and, hence, aperiodic electrohydrodynamic instability of the liquid 
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surface develops. If, however, the condition k > 1 holds, the frequency u is real-valued, 
which corresponds to the propagation of linear dispersive waves. 

The approach to the study of the evolution of a charged liquid surface based on the anal- 
ysis of relation (32) is obviously applicable only in the case of small-amplitude perturbations 
of the boundary: A <C k^ 1 . For finite-amplitude waves, the nonlinear effect may consist in 



the dependence of the dispersion relation on A (see, for example, 25]): 

uj = u>(k, A). 

The amplitude dependence of frequency is usually sought in the form of a power series in A 
(Stokes expansion), which limits the analysis to the weak-nonlinearity limit. Let us prove 
that, for electrocapillary waves, an exact solution to the nonlinear dispersion relation can 
be found. 

The equations describing a progressive wave (whose profile does not change in the refer- 
ence frame attached to the wave) can be obtained from the electrohydrodynamic equations 
(3)-(9) with the help of the following substitutions: 

(p — (f(x',z), $ = $'(x', z) + Cx', T] — T](x'), 

where x' = x — Ct and constant C has the meaning of the velocity of a wave moving in the 
direction of the This gives 

+ $' 2Z = 0, (33) 
(fx'x' + Vzz = 0, (34) 
2 + 2 "(1 + ^) 3 / 2 ' V{Xh (35) 

$; = ^$;„ z = v (x'), (36) 

V = 0, z = V (x% (37) 
$' -»■ -Cx', z -> -oo, (38) 
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y? — ► —z, z — > — oo. (39) 

These equations can be simplified by introducing the function of current ^(x',z), which is 
harmonically conjugate to potential 

This function satisfies the Laplace equation 

*W + *« = (40) 

with the boundary conditions 

$ = o, z = r?(x'), (41) 



* -> -Cz, z -> -oo, (42) 

which follow from relations (36) and (38). It can easily be seen that Eqs. (40)-(42) coincide 
with the Eqs. (34), (37) and (39) for the electric potential. Consequently, the following 
functional relation exists: 

* = Ca- 



using this relation, we can considerably simplify the Bernoulli equation (35), which assume 
the form 

C 2 + l ^ 2 , ,„ 2 ^ _ lh 

(1 + % 



— (^ + ^-1)= '/^ z = ^). (43) 



In combination with relations (34), (37), and (39), this condition completely defines the 
shape of a wave propagating in the coordinate system {x', z}. 

Equations (34), (37), (39) and (43) coincide except for constant factors with the equations 



describing the shape of a progressive capillary wave 



3 



and an equilibrium configuration of 



the charged surface of the liquid metal [26j. These equations have exact periodic solutions 
for which the boundary of the liquid is defined by the parametric expressions 

Ak~ 2 



2(C 2 + l) -1 + Acos(kp) 



+ z 0} (44) 



, 2Ak~ 1 sin(fcp) , 

X ~ P ~ 2(C 2 + l)- 1 + A cos(kp) + X °' (45j 
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where z and x' are constants, p is a parameter (the value of p changes over a period by 
27r/k), and the quantity A has the meaning of the amplitude of a surface perturbation; i.e., 
A = (zmax — £min)/2. The dependence of A on k and C is specified by the relation 

1/2 



.4 



(46) 



_(C 2 + 1) 2 fc 2 _ 

It was mentioned in Q| that solutions (44) and (45) exist only for 1 < k/(C 2 + 1) < 7, 
where 7 ~ 1.52. 

Considering that C is the phase velocity of the wave, we set C = ui/k in relation (46). 
Solving the obtained equation for frequency u), we arrive the exact nonlinear dispersion 
relation 

^ k ' A) =7mm- e ' (47) 

and the conditions of its applicability 

P7- 1 < u 2 - k 2 < k 3 . (48) 

It can be seen that, in the limit of infinitely small amplitudes (A —* 0), expression (47) is 
transformed into the linear dispersion relation (32). Let us consider the consequences of this 
nonlinearity. It can be seen from relation (47) that, for a fixed wave number k > 1, the 
maximum value of the surface perturbation amplitude A max (k) corresponds to the minimum 
possible value of u 2 . It follows from conditions (48) that, for 1 < k < 7, the value of 
^min = 0' which corresponds to a wave with zero velocity. In this case, expressions (44) and 
(45) define the solution of the problem on the steady-state profile of the charged surface of 
the liquid helium. For k > jk\, the amplitude has the maximum value for electrocapillary 
waves propagating at the velocity C = \fk^ x — 1; in this case, u>^ Din = k 3 ^ 1 — k 2 . This 
gives 

0, < k < 1 



A max (k) 



2y/l - k~ 2 , 1 < k < 7 
- 1, A;>7 

(the curve describing this dependence is presented in Fig. HJ). If the amplitude exceeds this 
value, expressions (44) and (45) describe a self-intersecting surface, which cannot be realized 
from the physical considerations, or w 2 < 0, which corresponds to incorrectly formulated 
problem in the context wave propagation. This leads to the assumption that the condition 
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FIG. 4: The maximum value of amplitude A max of an electrocapillary wave on the charged surface 
of liquid helium as a function of the wave number k. For k < 7, the peak corresponds to the value 
of uj = 0, while, for k > 7, the frequency differs from zero. 

v4(/c) > v4 max (fc) is the criterion of hard excitation of electrohydrodynamic instability of 
a plane charged suface of liquid helium, which generalizes the simplest linear instability 
criterion k < 1 to the case of finite-amplitude perturbation. 

It should be noted that peak of the function A max (k) corresponds to k — 7. The shape 
of the liquid surface corresponding to this value of the wave number is depicted in Fig. |3J 
It can be seen that the liquid acquires cavities. Such solutions reflect the tendency to 
the formation of charged bubbles (referred to as bubblons in the experimental work |12j |) 
on cuspidal dimples of the liquid helium boundary. The main mechanism of departure of 
electrons from the surface is associated with the generation of such bubbles. 
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FIG. 5: A period of the steady-state profile of a charged surface of liquid helium for k = 7. For 
this value of the wave number, the amplitude of the electrocapillary wave attains its maximum 
value. 

VI. AXISYMMETRIC SOLUTIONS 



Let us consider the evolution of the charged surface of liquid helium in an important case 
of the axial symmetry of the problem. The equations of motion (20)-(23) corresponding 
to the increasing branch of the solutions to system (10)— (15) in the cylindrical coordinaties 
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{r, z'} = {r,z — t} assume the form 

ip„. + r" Vr + </Vz' = 0, 

rj t = -(v? r 2 + + V r 2 )i = r/(r, 0, 

¥> = 0, z' = r)'(r,t), 

If — > — z', z' — > — oo. 

Here, r = x 2 + y 2 and we have taken into account the fact that <9 n y? = — |V</?| at the 
equipotential boundary in condition (21). 

At essentially nonlinear stages of the formation of a dimple on the surface of a liquid, we 
can assume that the electric field in the region of a large curvature of the surface is much 
stronger than the external field, |V</?| ^> 1. In this case, the dynamics of the boundary 
i]' = rj — t is completely determined by the intrinsic field rapidly attenuating with increasing 
distance, which allows us to use the condition |Vy?| — > for z — > — oo instead of the condition 
of field uniformity. We will also assume that the velocity of the liquid surface is considerably 
higher than the velocity of the origin in the laboratory reference frame (i.e., \i] t \ ^> 1). In 
this case, we can substitute 77 for r( and z for z' . This gives 

(p rr + r~ 1 (p r + (p zz = 0, z<r](r,t), (49) 

(p t = -(p r 2 - (p z 2 , z = r}(r, t), (50) 
V? = 0, z = rj(r,t), (51) 

ip r 2 + ip z 2 -> 0, r 2 + z 2 -> oo. (52) 
In relation (50), we have used the following conditions at the boundary of the liquid: 

Vt = -<Pt/<Pz, Vr = -¥r/Vz- 

The conditions of the applicability of this approximation will be considered at the end of 
this section. 
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A particular solution of Eqs. (49)-(52) can be obtain by using a substitution similar to 
that used in [2jj for constructing the axisymmetric solutions of the Stefan problem: 

V(r,z,t) = f(u(r,z,t)), (53) 



u(r, z, t)=-z-Vt+ yjr 2 + (z + Vt) 2 , (54) 

where the constant V has the meaning of the inward-directed velocity of the liquid surface. 
It can easily be seen that the equipotential surfaces corresponding to Eqs. (53) and (54) 
form a family of confocal paraboloids of revolution: 

r 2 = 2u(z + Vt) + u 2 (55) 

with the focus at the point r = and z = —Vt. 

Substituting expressions (53) and (54) into Eq. (49), we arrive at the following ordinary 
differential equation: 

ufuu + fu = 0. (56) 

It follows from Eqs. (50) and (51) that the boundary conditions to this equation have the 
form 

f u (u ) = V/2, f(u ) = 0. (57) 

Here, Uq is the value of parameter u at the surface of the liquid. Henceforth, we will use the 
quantity K = l/u which, in accordance with Eq. (55), defines the curvature of the liquid 
surface at the symmetry axis. Solving Eqs. (56) and (57), we obtain 

f(u)=Vln(Ku)/(2K), (58) 

which, together with Eqs. (53) and (54), describes the time evolution of the electric potential. 
It should be noted that condition (52) is naturally satisfied. The shape of the surface for 
the given exact solution of the equations of motion (49)-(52) is defined by the relation 

r](r, t) = Kr 2 /2 -Vt- (2K)-\ (59) 

which corresponds to needle-shaped dimple pulled into the bulk of the liquid velocity V. 
Such a geometry of the surface perturbation can be regarded as the simplest (paraboloidal) 
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approximation of the shape of the liquid boundary at essentially nonlinear stages of the 
development of instability of the charged boundary of the liquid. 

It should be recalled that the applicability of approximation (49)-(52) of the initial system 
(20)-(23) is limited by the conditions \r)t\ ^> 1 and |Vy?| ^> 1. Since rjt = —V in solutions 
(59) for any r and t, the first condition is reduced to the inequality V ^> 1 (in the dimensional 
notation, V ^> Ey/Anp). As regards the second condition, we can find the characteristic size 
D of the region in which the electric field created by a charged paraboloidal surface exceeds 
the external field. It follows from relations (53), (54) and (58) that the field distribution in 
the liquid is described by the relation 

' VV? ' ~ KV2Ru' 

Here, R = \Jr 2 + (z + Vt) 2 is the distance to the focus of the paraboloid; i.e., the field 
attains it maximum value equal to V at the point r = and z = —Vt — (2K)^ 1 . Since 
the field strength generally decreases in proportion to R^ 1 with increasing distance from 
the focus, the scale of D can be estimated as D ~ V/K. It should be noted that such a 
conclusion makes sense only if the value of D is much larger than the radius of the curvature 
K~ x of the surface perturbation. This requirement again leads us to the inequality V >> 1. 

Thus, we have obtained partial axisymmetric solutions to the equations of motion of 
liquid helium with a charged surface, which describe the evolution of a localized perturbation 
of the surface with a considerable curvature, and have determined the conditions of their 
applicability. However, the obtained solutions should not be regarded as general-position 
solutions. In all probability, solutions of the burst type, for which the surface becomes 
indefinitely cuspidate over a finite time interval, will dominate as in the 2D case. 

VII. CONCLUDING REMARKS 

In the absence of a surface charge, the equations of electrohydrodynamics of liquid helium 
considered by us are transformed into the well-known equations of a vortex-free flow of 
an incompressible liquid with a free boundary. These equations are extremely difficult to 
analyze, and the methods for the solution have not been developed at present. In this work, 
we succeeded in proving that the inclusion of the electrostatic pressure does not complicate 
the analysis of these equations. On the contrary, the emergence of an additional term in 
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the dynamic boundary condition introduces a certain symmetry into the equations so that 
they become compatible with the conditions <p + z = ±$. The emerging functional relation 
between the potentials of velocity and of electric field makes it possible to reduce by half the 
number of equations required for describing the motion of the surface and, in the long run, 
to find a wide class of exact solutions of the equations of motion of liquid helium with the 
boundary charged by electrons. It is important that the solutions obtained by us are not 
limited by the condition of smallness of surface perturbations; they describe the evolution 
of the liquid boudary up to the formation of cuspidal points in it. 

The dynamics of the formation of singularities in the case when the characteristic scale 
A of surface perturbations is comparable with the value of aE~ 2 and the capillary effects 
must be taken into consideration has not been considered by us here. In 2D geometry, such 
an analysis can be carried out using the methods of investigations of 2D potential flows 
with a free boundary, which was proposed in j3, and is based on conformal mapping 
of the region occupied by the liquid to a half-plane. In terms of the present work, such a 
transformation corresponds to the use of the field potential (p and its harmonically conjugate 
function v as independent variables. In the case of the axial symmetry of the problem (such 
a geometry reflects the experimentally observed phenomena most correctly), the for- 

mation of singularities can be described by self-similar solutions of the electrohydrodynamic 
equations, which are analogous to those considered in the recent publication |30| devoted to 
the formation of conic tips on the surface of a liquid metal in an external electric field. In 
accordance with the self-similar scenario of the development of instability, conical dimples 
with an angle of 98.6° appear on the surface over a finite time. A detailed analysis of these 
processes calls for further investigations. 

The author is grateful to V. E. Zakharov and E. A. Kuznetsov for their interest in this 
research. 
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